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Abstract 


We present an experimental study of the mixing proeesses in a gravity eurrent flowing on an inelined plane. 
The turbulent transport of momentum and density ean be deseribed in a very direet and eompaet form by 
a Prandtl mixing length model: the turbulent vertieal fluxes of momentum and density are found to seale 
quadratically with the vertical mean gradients of velocity and density. The scaling coefficient, the square of 
the mixing length, is approximately constant over the mixing zone of the stratified shear layer. We show how, 
in different flow configurations, this length can be related to the shear length of the flow We also 

study the fluctuations of the momentum and density turbulent fluxes, showing how they relate to mixing and to 
the entrainment/detrainment balance. We suggest a quantitative measure of local entrainment and detrainment 
derived from observed conditional correlations of density flux and density or vertical velocity fluctuations. 
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1. Introduction 

Mixing in stratified shear flows is an important 
process in many geophysical situations. Eerrari and 
Wunsch have detailed in 111] the kinetic energy bud¬ 
get occurring over the large range of oceanic scales. 
Several physical mechanisms, not all completely un¬ 
derstood, allow the energy to be transferred from 
large scale geostrophic motion to the very small 
scales where irreversible mixing can take place. Of 
particular current interest is the mixing and entrain¬ 
ment of oceanic overflows where gravity currents 
take place in particular regions of the oceans (Den¬ 
mark Strait flow, Mediterranean outflow). Although 
they occur in very localized regions, these currents, 
through their mixing processes, contribute strongly 
to the transport of heat and salinity in the riobal 
ocean via the thermohaline “conveyor belt” y, 

The mixing processes occur at scales too small to be 
captured by the numerical simulations of this circu¬ 
lation, requiring a sub-grid parametrization. In situ 
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measurements 0], as well as experimental studies, 
are necessary to provide an accurate description. In 
order to obtain a valid parametrization from a labo¬ 
ratory experiment, there is also a need for a model 
that extrapolates the parametrization to oceanic con¬ 
ditions. 

Ellison and Turner Isl] were the first to investigate 
this phenomenon experimentally, in the case of a cur¬ 
rent flowing into a homogeneous ambient medium, 
without rotation. They derived a model for the bulk 
properties of the flow, based on measurements over 
various quantities. Decades later, Baines 0,0] re¬ 
produced the experiment, with an ambient stratified 
flow, while more recently, Cenedese 0] studied grav¬ 
ity current in a rotating frame. 

In a broader context, several studies have been de¬ 
voted to the evolution of turbulence in a stratified 
shear flow. In such cases, the flow was generally 
not buOTancy-driven. Some of them, either numer¬ 
ically 0 or experimentally EQ, were focusing 
on general properties of turbulence, such as the in- 
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fluence of stratification and shear intensity on the 
growing or deeaying character of turbulence. Others 
looked at mixing properties in a similar way as we 
did HHQ. We give elsewhere ifTsIl a compari¬ 
son of our results with some of these studies. But in 
all these cases, for the experimental studies, the type 
of measurements used were either qualitative (video 
with dye) or quantitative but pointwise (conductiv¬ 
ity probe or density measured loeally on a fluid sam¬ 
ple). And in general, no veloeity measurement was 
performed. In order to better understand and model 
the proeesses involved in turbulent mixing indueed 
by gravity eurrents, it is essential to measure veloeity 
and density fields over a spatially extended region, as 
is now allowed by techniques such as particle image 
velocimetry and laser induced fluoreseence. 


We developed an experimental apparatus where an 
experimental gravity current ean be ereated with the 
ability to vary the flow eonditions sueh as the de¬ 
gree of stratification or the initial turbulence level. 
Among the key quantities for describing turbulent 
mixing processes are the fluxes of turbulent transport 
of momentum and density, which can be computed 
via the correlation terms between fluctuating compo¬ 
nents of the velocity field or between one component 
of velocity and density. One important feature of our 
experiment is the ability to measure simultaneously 
high resolution veloeity and density fields, thus al¬ 
lowing us to obtain sueh correlations. In the next 
section, we present the experimental set-up. Then, in 
section [3l we use our measurement of the turbulent 
fluxes to determine eddy diffusivities in our mixing 
processes. Section |4] recalls and extends the results 
of an earlier publication ifl^ . where we showed that 
a mixing length model provides a better description 
of how the turbulent fluxes are related to local mean 
gradients, compared to the constant eddy diffusivity 
hypothesis. In seetionO we use data taken in differ¬ 
ent parameter eonfigurations to determine a scaling 
law for the measured mixing lengths. Finally, in see- 
tion 0 we use the PDF’s of the fluxes to better un¬ 
derstand the location of entrainment and detrainment 
processes in the gravity current, before concluding in 
section |71 


2. Experimental set-up 


The experiment, sketehed in Fig. [H and described 
in detail elsewhere llisll . eonsists of a turbulent, 
uniform-density flow injected via a pump through a 5 
cm high by 45 cm wide nozzle into a tank filled with 
unstirred higher density fluid. The turbulence level 
of the injection current can be enhanced by an ac¬ 
tive grid device located just before the injection noz¬ 
zle. The flow, upon exiting the nozzle, is bounded 
from above by a transparent plate inclined at an an¬ 
gle of 10° with respeet to horizontal, is unbounded 
below, and is confined in a tank about 2 m long, 
0.5 m wide and 0.5 m high. The components of the 
spatial position vector x describing the flow are the 
mean flow direction x, the cross-stream direetion y 
and the downward distance perpendicular to the plate 
z. The corresponding velocity u(x) has components 
{u,v,w}. We use the notation (u) for a time- and 
ensemble-averaged quantity and u' = u - (u) for its 
fluctuating paro The injection fluid, a solution of 
ethanol and water, is less dense than the fluid in the 
tank, water and salt (NaCl). This situation is reversed 
compared to oeeanic overflows, where a denser fluid 
flows down an ineline, but within the Boussinesq 
approximation, the physieal mechanisms involved 
are the same. The density differenee is defined as 
Ap = p - ps, where p is the measured density and p^ 
is the maximum density, corresponding to the initial 
salt water density. The eoncentrations of ethanol and 
salt are adjusted (and the fluid temperatures main¬ 
tained equal within 0.2°C) so that the fluids are index 
matehed to avoid optieal distortions IllSll . Eaeh eon- 
figuration presented in this study consisted in several 
runs with the same eonditions, adding up to a total 
between 500 and 800 veloeity fields. All the fluids 
are freshly prepared for eaeh run. 

Instantaneous veloeity and density fields are mea¬ 
sured in a 9 cm x 9 cm area of a 0.1 cm thick 
laser sheet in the x - z plane. Velocity and den¬ 
sity are measured simultaneously using particle im¬ 
age velocimetry (PIV) and planar laser-induced fluo¬ 
rescence (PLIF), respectively, at a rate of 3 Hz with 


'in our previous publication, 13], ■ represents ensemble 
averaging and (■) denotes spatial averaging. This adjustment in 
nomenclature is made to be consistent with the popularly used 
ones in literatures, e.g. IS. 
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two 2048^ pixels digital cameras. Fluorescent dye 
(Rhodamine 6G) is added to the light fluid, and a 
calibration of density versus fluorescence intensity is 
performed for each position of the field of view. 

The lighter exit fluid is stably stratified with re¬ 
spect to the heavy fluid in the tank and forms a grav¬ 
ity current on the bottom side of the plate. The 
competition between the stabilizing effect of buoy¬ 
ancy and the destabilizing shear is captured in a 
dimensionless parameter, the Richardson number, 
Ri = -(gAplpsH)l(U^) where g is the acceleration 
of gravity, U a typical velocity, and H the height 
of the injection nozzle. For small Ri, shear domi¬ 
nates buoyancy, and the flow is unstable to Kelvin- 
Helmholtz instability lll9n . In a standard configura¬ 
tion, used for the results in section [3l |4] and the 
current is injected at a speed of Uq = 1 cm/s, with the 
active grids on, and the initial density difference be¬ 
tween the fluids is Apo = 2.6 g/L. This results in an 
initial value of Ri = 0.27 and a fully turbulent flow as 
it exits the nozzle with streamwise velocity fluctua¬ 
tions u' about 25% of (w), corresponding to a Taylor 
Reynolds number Re^ = u'^! VlSev « 100, where 
V is the fluid kinematic viscosity and s is the mean 
dissipation rate measured directly from velocity field 
(the spatial resolution of our velocity measurement 
is 0.5 mm compared to the dissipation scale of 0.33 
mm). 

Other configurations (see section [S]) were tested, 
where injection speed and/or initial density differ¬ 
ence were varied, and in one case the active grids 
were removed. 

In all configurations, there is rapid evolution of 
mean quantities over the first 20 cm. In this report, 
we focus on the region from 21 to 45 cm over which 
averages are approximately uniform along x. Note, 
however, that the results described here also apply 
to the initial region, except with a stronger depen¬ 
dence on downstream distance. A detailed descrip¬ 
tion on the downstream evolution of various quanti¬ 
ties associated to the current has been reported else¬ 
where lll5n . 



Figure 1: Sketch of the experimental device. 

measure in our experiment, can be interpreted, re¬ 
spectively, as the vertical^ flux of downstream mo¬ 
mentum and of density due to turbulent fluctuations. 
We show in Fig. |2] the vertical profiles of these two 
quantities computed locally as an ensemble average 
over all PIV/PLIF images and over all experimental 
runs. As expected, they display a maximum in the 
mixing region, close to the initial interface between 
the current and the ambient fluid. Close to the plate 
the mixing is reduced because perturbations that ad- 
vect high concentration regions towards the plate are 
very rare, and far from the plate the ambient fluid is 
also undisturbed. In addition, the amplitude of the 
fluxes decreases as one goes away from the injection 
nozzle. This reduction reaches a factor of 1.5 for the 
momentum flux when the distance is doubled, and a 
factor of 3 for the density flux. 

In models of oceanic circulation, various closure 
schemes are used to parametrize turbulent transport 
(for a review, see [^, ^). Some of these models are 
highly elaborate, closing the turbulence equations ei¬ 
ther at first moment level, providing sets of differ¬ 
ential equations relating turbulent fluxes to the mean 
motion or even at second moment level, with equa¬ 
tions for quantities like pressure strain-rate correla¬ 
tions. In more basic models (zero-equation models), 
the relation between turbulent fluxes and mean gra¬ 
dients is considered to be an algebraic one. For ex¬ 
ample, the model initially proposed by Pacanowski 
and Philander ll22n . assumes a proportionality be¬ 
tween fluxes and gradients, the proportionality con¬ 
stant being an effective diffusivity times a simple an¬ 
alytical function of the Richardson number. This 
proportionality constant is thus defined as vj = 
-{u'w') I(d^u) in the case of momentum transport and 
7t = -{p'w')l{d-p) for density transport. 


3. Eddy diffusivities 

As mentioned in the introduction, correlation 
products such as {u'w') and {p'w'), which we can 


^Note that for simplicity we use the word “vertical” for a 
direction that, strictly speaking, makes a 10° angle with the ver¬ 
tical. 
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Figure 2: Vertical profiles of the measured momentum turbu¬ 
lent flux (a) and density turbulent flux (b) measured at differ¬ 
ent distances downstream. 

Our respective measurements of the fluxes and the 
corresponding gradients allow us to make an exper¬ 
imental determination of the eddy diffusion coeffi¬ 
cients. Figure [3] shows the vertical profiles of the 
measured eddy diffusivities, in the case of momen¬ 
tum transport (coefficient vt) and density transport 
(coefficient jt). They are not shown for extreme lo¬ 
cations, too close to the plate (x < 1.5 cm) and too far 
away from the plate (x > 7.5 cm), where the mixing 
is too weak to allow reliable measurement.They are 
displayed for 3 different distances from the injection 
nozzle. 

It is interesting to note that the measured value 
of both coefficients is the same, about 30 times the 
molecular viscosity. In addition, these values vary 
by a factor of 3 between the central region, where 
most of the mixing occurs, and the top and bottom 
regions, which are respectively the less disturbed re¬ 
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Figure 3: Vertical profiles of the measured eddy viscosity 
(a) and eddy density diffusivity (b) computed at different dis¬ 
tances downstream. 


gion in the current near the inclined plate and the less 
disturbed region in the ambient fluid away from the 
current. The coefficients also seem to decrease as 
one measures them further away from the injection 
nozzle. 

For comparison, we show in Fig. |4] the corre¬ 
sponding plot (analogous to Fig. [3^) in a case where 
there is no density difference between the injected 
fluid and the ambient fluid. This configuration cor¬ 
responds to a plane jet along a wall. One can ob¬ 
serve that the dependence of the eddy viscosity with 
the distance to the wall is much weaker: the verti¬ 
cal variation is of the order of 30%, instead of the 
factor 3 observed in the stratified case. This observa¬ 
tion is consistent with what has been measured in the 
case of unstratified jets (see for example figure 5.10 
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in H). In the same referenee, the study of a plane 
self-similar jet shows that the eddy viscosity is ex¬ 
pected to increase with downstream distance, which 
is what we observe in Fig. HI but not in Fig. [3l That 
the constant eddy viscosity model works better in the 
unstratified flow and fails rather badly in our strati¬ 
fied flow experiments is interesting but we are unable 
to explain this difference. 


small enough so that the variation of the gradient 
over a distance can be neglected, one can take 
f/(typical) = L,„\{d^u)\ and thus obtain the relation 
{u'w') = LlAdyU)^, which is indeed observed (see 
figure 3 in niblF . The same argument for the density 
flux yields: {p'w') = -Lp\{d^u)\{d;p), where Lp is a 
mixing length associated with the density transport. 

As a result, we computed the mixing lengths as: 
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Figure 4: Vertical profiles of the measured eddy viscosity, 
computed at different distances downstream, for an unstrat¬ 
ified flow. Because of the very different downstream evolu¬ 
tion in this case, some measurements are shown for distances 
downstream different from the ones used in the other figures, 
and are therefore shown with different colors and symbols. 


4. Mixing length model 

The strong dependence of the turbulent diffusivi- 
ties with depth and distance downstream calls for a 
better scheme to model the relation between fluxes 


Ll = 


{u'w') 

{d,u)^ 


and Lp = 


-{p'w') 


mu)\{dj)) 


(1) 


and mean gradients. In an earlier publication [Il6n . 
we showed that the observed scaling between these 
quantities is actually : {u'w') oc {d^u)^ and {p'w') oc 
\{d^u)\d^. We also demonstrated that this observa¬ 
tion can be understood in the framework of Prandtl 
mixing length theory ll^ . Prandtl’s argument is 
analogous to that applied in the kinetic theory of 
gases to molecular transport processes: it assumes 
that the coefficient of eddy viscosity is equal to 
the product of a “mixing length” L,„, characteris¬ 
tic of the mixing phenomena, by a suitable veloc¬ 
ity: vt - X f/(typical). Assuming that L,„ is 


The resulting vertical profiles of mixing lengths 
are shown in Fig. |5l They are much more uniform 
over depth, compared to the turbulent diffusivities 
profiles shown in Fig. HI The dependence in down¬ 
stream distance is also much weaker. No decrease of 
the mixing lengths is observed as one approaches the 
wall (z = 0), but one must remember that a law of the 
wall should only start within a distance from the wall 
of approximately 20% of the flow width [Il7n . about 
1 cm in our case, which is a region where we have no 
accurate measurement of the mixing lengths (see the 
remark about measurement errors in the caption of 
Fig. O. The mean values (average taken over depth 
and downstream distance) for both mixing lengths 


are the same: L„ 


5. Shear scale 


0.45 + 0.1 cm 


In order to allow extrapolations of these results to 
oceanic conditions, it is necessary to understand how 
the mixing lengths depend on the degree of stratifi¬ 
cation of the current, as well as on its level of tur¬ 
bulence. In lll6n . we derive an interpretation of the 
mixing length using the balance between the produc¬ 
tion of turbulent kinetic energy by the destabilizing 
shear and its dissipation. This balance defines a scale 

Ls = (sKd^u) which we call “shear scal e” 11241] . 
This scale was defined initially by Corrsin I125n as 
the smallest scale at which anisotropy effects result¬ 
ing from a large scale shear are carried out by the 
turbulent cascade. This scaling has been confirmed 
experimentally in a boundary layer flow ll^ . As 
mentioned earlier, at low enough Richardson num¬ 
ber, the effect of shear dominates the effect of buoy¬ 
ancy, therefore the relevant quantity to define the 
mixing scale is the shear and not the Brunt-Vaisala 
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Figure 5: Vertical profiles of the measured mixing lengths for 
momentum (a) and density (b) at different distances down¬ 
stream. Note that close to the plate, as well as far away from it, 
the gradients tend to vanish, producing large errors in the mix¬ 
ing length calculation. This explains the “stray” data points 
which can be observed at the extremity of some curves. 


frequency N = Since we are inter¬ 

ested in a global scale for the flow, a spatial aver¬ 
age (“=“ symbol) is taken for the quantities used to 
compute this scale, which corresponds to a different 
definition than the one used for the mixing lengths 
in equation [TJ This is important because the whole 
point of this section is to show that since the mixing 
lengths, although they are defined locally, do not de¬ 
pend on the local position, they can be related to a 
global quantity, namely the shear scale. 

We show in [Q that the value of the measured 
mixing lengths in the standard configuration is very 
close to the shear scale. In order to test the robustness 


of this observation, we measured the mixing lengths 
in various flow configurations. The results are dis¬ 
played in table [B Configuration 1 (shaded line in the 
table) is the standard configuration used until now 
(and also used for the data presented in lll6l] l. We 
then vary stratification and/or intensity of turbulence 
in 5 different configurations. Since in all cases mo¬ 
mentum and density mixing lengths are equal (ex¬ 
cept the non stratified case where Lp is undefined), 
we give in column 7 the value of L,„ only. 

Compared to case 1, cases 2 and 3 show that the 
value of the mixing length increases with turbulence 
intensity^ In the same way, as expected, cases 4 and 
5 show that the stratification prevents mixing. In the 
last case, both the stratification is stronger and the 
shear weaker, resulting in a very short mixing length. 

In Fig.O we plot the measured mixing length ver¬ 
sus the computed shear scale, showing that in all 
cases studied, there is a scaling between the two 
quantities: L„, = (0.7 + 0.03)Ls. In the unstratified 
case, we observed that the quadratic scaling between 
momentum flux and velocity gradient is less evident, 
which may be related to the observation that the con¬ 
stant eddy viscosity model seems to work better in 
this case, as shown in section [3l We included this 
data point in the graph however, since interestingly 
enough, it fits nicely the general trend, with a much 
larger value for the scales. One can see from the in¬ 
set, which shows an expansion of small lengths (cor¬ 
responding fit indicated in the main plot by a dashed 
line), that the slope is not very different (5% change) 
when the unstratified point is included (or not) in the 
fit. 

It may seem surprising that the mixing length 
scales with a quantity apparently independent of the 
stratification. But one must keep in mind that the 
mean turbulent dissipation rate e, and the mean shear 
(d^u) depend on the stratification. A study we per¬ 
formed, presented elsewhere lllSn . shows that the 
shear is increased by the presence of stratification, 
because this stratification partially prevents the mix¬ 
ing that reduces the shear. This observation is in 


^Note, however, that the quadratic relation between momen¬ 
tum flux and velocity gradient could not be clearly observed in 
the unstratifled case because of a much larger spread of the data 
points, possibly owing to stronger mixing. This makes the cor¬ 
responding calculation of the mixing length less reliable. 
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# 

configuration 

Apo/p.9 

Ri 

t/o 

[cm/s] 

Ra 

Lfn 

[cm] 

(s) 

[cm^/s^] 

[s-'] 

Ls 

[cm] 

1 

standard 

0.26% 

0.3 

7 

100 

0.45 

0.8 

1.4 

0.55 

2 

larger injection speed 

0.26% 

0.2 

9.5 

140 

0.6 

1.4 

1.4 

0.7 

3 

no active grid 

0.26% 

0.35 

6.5 

42 

0.35 

0.7 

1.7 

0.35 

4 

unstratified 

0 

0 

7.5 

115 

2.1 

1.5 

0.6 

2.7 

5 

double density 

0.52% 

0.45 

7 

93 

0.3 

0.8 

2.1 

0.3 

6 

double dens, 
half veloc. 

0.52% 

1.4 

4.3 

72 

0.2 

0.5 

2.7 

0.15 


Table 1: Summary of the different experimental configurations. The first column is a reference number, the second gives the general 
features of the configuration, compared to the standard one (shaded line). The third column gives the initial density difference and 
the fourth the Ri. The fifth column shows the initial velocity of the current and the sixth shows R,i. The 3 last columns give the data 
necessary to compute the shear length, and the shear length itself. The mean shear (column 9) is averaged over the half width of 
the vertical profile of the gradient. 


agreement with the observed deerease of the mixing 
length when the stratification is stronger. In the case 
of the turbulent dissipation rate, its dependence on 
the stratification is more subtle, since the buoyancy 
force acts at the same time as a turbulence suppres¬ 
sor and as a shear enhancer, since it is the driving 
force of the current. 



Figure 6: Measured mixing length vs computed shear length, 
for the different flow configurations indicated. The error bars 
are determined by the standard deviation of the mixing length 
values at different locations in the flow. The solid line is a linear 
fit to the data. The inset shows a expanded view of the region 
below 1 cm, with a fit not including the point for the unstratified 
case. The slope for this fit is also shown in the main plot by a 
dashed line. 


6. Fluctuations of turbulent fluxes and entrain¬ 
ment 

The mixing length model and the scaling proper¬ 
ties presented in the previous sections allowed a bet¬ 
ter understanding of the average behavior of the tur¬ 
bulent momentum and density fluxes. In this sec¬ 
tion, we focus on the fluctuations of these fluxes, 
and how they relate to the mechanisms of entrain¬ 
ment and detrainment of the flow. Entrainment oc¬ 
curs when a parcel of fluid that is heavier (in our 
case) than the fluid within the gravity current and lo¬ 
cated outside that current is advected into the cur¬ 
rent and thoroughly mixed. Average measures of 
entrainment show how much the gravity current in¬ 
creases in overall volume flow rate. Similarly, de¬ 
trainment is the process in which a lighter element of 
fluid is transported by turbulent fluctuations into the 
heavier fluid outside the gravity current and is ab¬ 
sorbed by mixing. Below we present a novel method 
for characterizing the distribution of correlations of 
density and vertical velocity fluctuations that con¬ 
tribute to understanding the dynamics of the entrain- 
ment/detrainment process. 

The probability density function (PDF) of the mo¬ 
mentum and density fluxes, {u'w') and {p'w'), re¬ 
spectively, are shown in Fig. |71 We have to remind 
the reader that using our convention of signs for w 
and p, positive w' corresponds to a downward veloc¬ 
ity fluctuation and positive p' corresponds to a fluctu¬ 
ation of a particle heavier that the local mean. Each 
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color corresponds to a given horizontal band, thus 
allowing one to see the evolution of the PDFs as z 
inereases. First, the fluxes reaeh large values eom- 
pared to the mean as shown in Fig. [2l there is still a 
probability 10~^ that a fluctuation will reach a value 
about 10 times the mean. As expected, this tendency 
to produce large fluctuations is stronger as one ap- 
proaehes the mixing region around the initial inter- 
faee between the eurrent and the ambient fluid (blaek 
eurve, z=3-4 em). 

The PDFs of momentum and density fluxes are 
asymmetrie. This asymmetry is the origin of the non¬ 
zero mean value of the fluxes. With our eonventions 
of sign, it is positive for momentum and negative for 
density. In order to understand this asymmetry, it is 
neeessary to give some eonsiderations to the signs of 
the fluxes. As ean be seen in Fig.[Tl the x axis is ori¬ 
ented in the direction of the flow and the z axis points 
downwards. Thus, a positive value of u'w' eorre- 
sponds to either downward transport of downstream 
momentum or upward transport of upstream momen¬ 
tum. It is therefore understandable that the PDF of 
u'w' displays more positive events, sinee these eorre- 
spond to the standard transport of momentum owing 
to Kelvin-Helmholtz mixing. In other words, down¬ 
stream fluetuations with downward transport take 
turbulent momentum from where it is large, namely 
in the interior of the mixing zone, towards the qui- 
eseent regions in the heavier fluid below the gravity 
eurrent, i.e., at large z. In the same manner, negative 
values of p'w' eorrespond either to upward transport 
of heavier fluid or to downward transport of lighter 
fluid. Either of these transport meehanisms oppose 
buoyancy and are caused again by Kelvin-Helmholtz 
mixing. This explains the negative asymmetry ob¬ 
served in Fig. |7j5. In addition, we observe that this 
asymmetry becomes stronger for the PDF in the een- 
ter region (red and green, then black curve), where 
most of the mixing takes place. Finally, the PDFs 
of both momentum and density are broadest for the 
blaek eurves, i.e., in the middle of the mixing zone. 
On the other hand, the most extreme events oeeur on 
either side of the peak of the mixing zone where the 
average fluxes are largest. 

In order to eompare the momentum and density 
flux PDFs to more standard probability distributions 
studied in turbulenee, we show in Fig. [8] PDFs of 




Figure 7; PDFs of the momentum flux (a) and density flux (b). 
Each PDF is constructed using data in a horizontal band of 1 
cm height (vertical position indicated in the legend), situated 
between 20 and 49 cm from the injection nozzle. 

the vertieal veloeity fluetuations, normalized by their 
rms value, for horizontal layers at different distanees 
from the plate. The PDFs are fairly Gaussian in the 
mixing region (blaek line, eompared to a gaussian 
fit, indieated by a dashed line), whereas elose to the 
plate or far away from it, the distributions differ from 
a Gaussian shape, with broader tails, probably ow¬ 
ing to the intermittent and not fully-developed turbu¬ 
lence in these regions. 

One ean further elueidate the nature of the den¬ 
sity flux by eonsidering the instantaneous eorrela- 
tions between density fluctuations p' and density flux 
p'w'. To understand these eorrelations, one needs to 
eonsider the expected behavior of a pareel of fluid 
in a baekground stratification defined by the average 
density p(z). In the absenee of turbulenee, a pareel 
of fluid with positive or negative p' should move to- 
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Figure 8: Normalized PDFs of the vertical velocity fluctua¬ 
tions. Each PDF is constructed using data in a horizontal band 
of 1 cm height (vertical position indicated in the legend), sit¬ 
uated between 20 and 49 cm from the injection nozzle. 


wards its neutral buoyancy position, exchanging its 
potential energy for kinetic energy. In the case we 
consider, a parcel of fluid that is lighter than its sur¬ 
roundings has a negative value of p' and should move 
towards the plate, that is, have a negative vertical ve¬ 
locity w' < 0 so that the resultant density flux is pos¬ 
itive p'w' > 0, whereas a parcel of heavier fluid has a 
positive p' and would be expected to have a positive 
vertical velocity so again one has positive density 
flux, p'w' > 0. This situation is illustrated schemat¬ 
ically in Fig. and indicates that for either positive 
or negative p', stabilizing return to neutral buoyancy 
corresponds to positive density flux. Here, and in fur¬ 
ther discussion below, we use an intuitive description 
based on tracking a fluid parcel in time, i.e., a La- 
grangian perspective. Our measurements, however, 
are Eulerian so that the parcel of fluid we consider 
at time t is advected away by the mean flow. In ad¬ 
dition, the buoyancy condition implied by specifying 
a value of p' is not the instantaneous buoyancy ex¬ 
perienced by a fluid parcel because p' is determined 
by a long-time density average rather than by the in¬ 
stantaneous distribution of density. Nevertheless, the 
description presented here gives a reasonable under¬ 
standing of our results and provides the basis for fur¬ 
ther exploration of these important correlations. 

If we add turbulent fluctuations, the situation is a 
bit more complex, see Fig.l^^. A positive vertical ve- 


a) Stabilizing 

w' < 0 


p 


P' p'w'>0 
w' > 0 


b) Destabilizing 


w' > 0 
p'w' < 0 ^ 

^ p' <0 


P' > 0 , , ^ 

X- - kaP w'<0 

w' <0 



Figure 9: Schematic illustration of a) stabilizing velocity and 
density correlations resulting in p'w' > 0, b) destabilizing ve¬ 
locity and density correlations with p'w' > 0, and c) the net 
effect of mixing on fluid parcels that mix during a destabilizing 
displacement before gravitational forces can restore the parcel 
to its equilibrium vertical position of neutral buoyancy. In c), 
the solid (dashed) arrows in c) indicate non-mixing (mixing) 
processes, and the solid (dashed) curves show the non-mixing 
(mixing) mean density profile. 
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locity fluctuation w' > 0 may advect a lighter parcel 
of fluid with p' < 0 into a region of higher density, 
thereby producing a negative density flux p'w' < 0 
and a condition which is unstable with respect to the 
potential energy of that parcel. Similarly, a negative 
value of w' < 0 advecting a parcel of fluid with p' > 0 
also produces an unstable condition with p'w' < 0. 
Thus, destabilizing fluctuations have negative den¬ 
sity flux p'w' < 0. 

If no mixing occurs, these two processes of desta¬ 
bilizing and stabilizing motions should balance in 
that a parcel that is displaced vertically by turbulent 
fluctuations would be restored to neutral buoyancy 
owing to its adverse potential energy with respect to 
the mean density profile. This process is illustrated in 
Fig. [9}: by the solid arrows. If mixing occurs along 
this path, however, there will be a net effect on the 
local mean density profile that reduces the local den¬ 
sity gradient as indicated in Fig. by the dashed 
arrows; a positive vertical velocity (downward) pro¬ 
duces a net reduction in density upon mixing, i.e., 
6p < 0 whereas an initially negative w will produce a 
net increase in density, 6p > 0. The overall effect of 
many such parcels being displaced away from neu¬ 
tral buoyancy, mixing partially or completely, and 
returning to a new neutral buoyancy vertical position 
is to reduce the global mean density gradient. These 
examples described in Fig. ^ allow a reasonable in¬ 
terpretation of the experimentally obtained correla¬ 
tion PDFs between density fluctuations and density 
flux shown in Fig. (TOl In particular, it allows one to 
obtain the amount of entrainment (fluid added to the 
gravity current) or detrainment (fluid left behind in 
the ambient fluid) of the flow as a function of verti¬ 
cal height. The balance between entrainment and de¬ 
trainment in a boundary flow is an unresolved issue 
in numerical simulations and may be resolved 
using the correlations between p' and p'w'. 

Fig.fTOlshows 2D PDFs of p' versus p'w', in each 1 
cm horizontal band taken at different distances from 
the inclined plate. From the earlier discussion illus¬ 
trated in Fig. |9l we can tell that the right part of the 
plots in Fig. [TOl corresponds to the stabilizing return 
to neutral buoyancy (positive flux, see also Fig. |9^), 
while the top left quadrant corresponds to entrain¬ 
ment (see right part of Fig. [^5) and the bottom left 
quadrant to detrainment (see left part of Fig. I^^). 


One observes that in the regions close to the plate, the 
distributions are fairly symmetric with respect to the 
sign of the fluxes, with only a slight negative skew¬ 
ness in the positive p' lobe, indicating a small level 
of entrainment. Moving away from the plate, the 
asymmetry in the upper entrainment lobe increases, 
i.e., more entrainment, and the detrainment lobe with 
p' < 0 begins to become asymmetric as well. Near 
the middle of the mixing zone with z ~ 3 - 4 cm, 
the entrainment and detrainment lobes are about the 
same, indicating an approximately equal amount of 
entrainment and detrainment. Finally, far from the 
plate entrainment is inactive as indicated by the small 
and symmetric entrainment lobe compared to ac¬ 
tive detrainment demonstrated by the large and nega¬ 
tively skewed detrainment lobe. The qualitative pic¬ 
ture one gets from this analysis is that deep in the 
gravity current entrainment dominates whereas out¬ 
side the current detrainment is dominant. 

The correlations presented in Fig. [TO] should also 
be reflected in correlations between vertical velocity 
fluctuations and density flux. In Fig. [TTl we show a 
comparison between p'w'\p' correlations in Fig. [TOh 
and p'w'\w' correlations for conditions near the plate, 
z = 0 - 1 cm. One gets a similar skewness indicat¬ 
ing entrainment (the quadrants for entrainment and 
detrainment are flipped vertically for w' compared 
to p') although there are long tails in the stable re¬ 
gion of the p'w'\w' correlations. The precise details 
of the shapes of these PDFs is beyond the scope of 
this paper but overall they point to an interesting ap¬ 
proach at a local measure of entrainment and detrain¬ 
ment. In particular, one can compute, based on the 
analysis above, an effective entrainment parameter 
Ce = -{p'w'\p' > 0) and a corresponding detrain¬ 
ment parameter Cu = -{p'w'\p' < 0) (minus sign 
added because the fluxes are negative). In Fig. 
we show the quantities Ce and Co as functions of 
z. The results show the behavior described above, 
namely that entrainment is dominant near the plate 
whereas detrainment is larger in the quiescent fluid 
far from the plate. The combined effect of entrain¬ 
ment and detrainment is described by the addition of 
the two parameters, i.e,. Ce + Co, shown in Fig. O 
as (Ce + Cd)I2. The maximum density flux happens 
in the middle of the mixing zone, between 3 and 5 
cm. 
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Figure 10: 2D PDFs of the density flux vs density fluctuations. Each PDF is constructed using data in a horizontal band of 1 
cm height (vertical position indicated in the top right corner of each plot), situated between 20 and 49 cm from the injection 
nozzle. 11 
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Figure 11: 2D PDFs of a) density flux vs density fluctuations 
and b) density flux vs vertical velocity fluctuations. Each PDF 
is constructed using data in a horizontal band of 1 cm height 
next to the inclined plate, in a region between 20 and 49 cm 
from the injection nozzle. 
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Figure 12: Entrainment and detrainment parameters Ce (cir¬ 
cles), Co (squares), and their average value {Ce + Co)/2 (trian¬ 
gles) as functions of z. 


7. Conclusion 


In order to allow a comparison of our results 
with real oceanic situations, some additional features 
would need to be taken into account. To start with, 
as in any experimental study of an oceanic process, 
the Reynolds number dependence needs to be inves¬ 
tigated in more detail, since this number in the ocean 
is far greater than can ever be achieved experimen¬ 
tally. In addition, although the scales at which the 
oceanic mixing takes place are much smaller than the 
Rossby deformation radius, it has been shown [l28ll 
that the combined effect of the Coriolis force and 
the presence of topography (ridges or canyons on the 
bottom surface) can lead to secondary flows that can 
contribute to the mixing. Even without the Coriolis 
force, secondary flows can be created by topographic 
effects, for example in a curved channel. Finally, the 
roughness of the bottom (or of the inclined plate in 
our case) has been shown to influence the behavior of 
gravity currents [|3]- Extensions of this study, taking 
some of these aspects into account, would contribute 
to possible extrapolations to oceanic situations. Cer¬ 
tain conclusions can already be drawn, however, re¬ 
garding mixing processes in gravity currents. 

Our experimental device provides us with the abil¬ 
ity to derive correlations between components of 
velocity, as well as between velocity and density. 


12 














thanks to simultaneous measurements of these quan¬ 
tities. In addition, good spatial resolution allows us 
to study the evolution of these fluxes along both the 
flow direction and the vertical direction perpendic¬ 
ular to the plate. We measure eddy viscosity and 
diffusivity, although the constant eddy diffusivity as¬ 
sumption does not work well for the data. A mixing 
length model provides a better fit of the data with the 
definition of a typical length associated with the mix¬ 
ing phenomena. Contrary to the eddy diffusivity, this 
mixing length is very constant in space, allowing us 
to study its scaling with a length scale constructed 
from the turbulent dissipation rate s and the mean 

_ 3 

shear, = (sKd^u) Using data taken in differ¬ 
ent configurations of turbulence level and stratifica¬ 
tion, we show that this scaling is robust, even in the 
case of a simple wall jet (no stratification), where the 
mixing is much stronger with a mixing length almost 
ten times larger than in the stratified cases. 

Finally, mixing events can also be observed by 
looking at the asymmetry of the probability density 
functions of the momentum and density fluxes. As 
expected, the strongest mixing takes place at the in¬ 
terface between the current and the ambient fluid. In 
addition, 2D PDFs of the correlation between the 
density flux and the density fluctuations provide a 
local measure of entrainment and detrainment. In 
particular, we demonstrate that close to the plate 
the mixing is predominately entrainment of heav¬ 
ier fluid, whereas away from the plate, the mixing 
is largely associated with detrainment, where the 
current releases some of its fluid into the ambient 
medium. Quantitative determinations of entrainment 
and detrainment are derived from the asymmetry of 
the PDFs and provide a possible means for obtaining 
a global measure of net gravity current entrainment 
from local measurements. Connecting these results 
with more traditional measures of entrainment, i.e., 
mean vertical velocity divided by mean downstream 
velocity would be an very interesting extension to the 
present work. 
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Aluie, J.-F. Pinton and M. Rivera. Work performed 
at Los Alamos National Laboratory was funded by 
the US Department of Energy under Contract No. 
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